Evaporation of buffer gas-thermalized anions out of a multipole rf ion trap 
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We identify plain evaporation of ions as the fundamental loss mechanism out of a multipole ion 
trap. Using thermalized negative Cl~ ions we find that the evaporative loss rate is proportional to a 
Boltzmann factor. This thermodynamic description sheds new light on the dynamics of particles in 
time-varying confining potentials. It specifically allows us to extract the effective depth of the ion 
trap as the activation energy for evaporation. As a function of the rf amplitude we find two distinct 
regimes related to the stability of motion of the trapped ions. For low amplitudes the entire trap 
allows for stable motion and the trap depth increases with the rf field. For larger rf amplitudes, 
however, rapid energy transfer from the field to the ion motion can occur at large trap radii, which 
leads to a reduction of the effective trapping volume. In this regime the trap depth decreases again 
with increasing rf amplitude. We give an analytical parameterization of the trap depth for various 
multipole traps that allows predictions of the most favorable trapping conditions. 

PACS numbers: 05.20.-y,32.80.Pj,41.90.+e 
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Evaporation of atoms and molecules out of a confined 
thermalized ensemble is a well understood process [l|, [1] 
and represents the decisive cooling step towards Bose- 
Einstein condensation [1, In radiofrequency (rf) ion 
traps, early studies of the dynamics have shown that, in 
the absense of a buffer gas, thermalization is acchieved 
by ion- ion collisions and evaporative losses Q . In ap- 
plications of rf Paul traps [3| for quantum information 
processing [1, [9| , precision spectroscopy [l^, [ll|, , and 
the pro duction of translationally ultracold molecular ions 
[TS. Il4|. trapped ions are cooled using light forces. For 
the sympathetic cooling of ions in a buffer gas, higher 
order multipole rf traps at cryogenic temperature have 
proven to be more useful due to their large field free 
region [l^, [l^. Multipole ion traps are hence widely 
used to prepare cold molecular ions in laboratory astro- 
physics [l3] , to reduce the Doppler shift in microwave ion 
clocks '18|, to measure absolu te p hotodetachment cross 
sections of molecular anions [19|, the photofragmenta- 
tion of biomolecules (20| , for precision rovibrational spec- 
troscopy Hi], [23, [2^ and for collision experiments [Ij . 

If the ion temperature is fixed by collisions with a 
buffer gas, evaporative loss measurements provides access 
to the stability of ion motion in traps. In a quadrupole 
or Paul trap, stable ion motion is described analytically 
by the Matthieu equations, which puts well-defined con- 
straints on the trapping fields. In addition to these global 
boundaries of stability, heating of trapped ions by the rf 
field might be assisted or damped by Coulomb interac- 
tion, and phase transitions between crystalline and chaot- 
ically moving clouds are observed [25j . In contrast, for 
ions moving in a high order multipole field the equations 
of motion have no analytical solution. There have been 
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attempts to numerically establish a stability diagram for 
two-dimensional hexapoles and octupoles , 26;] . For differ- 
ent starting conditions the obtained stability diagrams 
look very different and the regions of stability do not 
show well defined boundaries opposite to the case of the 
quadrupole trap. 

In this letter we show that one can attribute in a ther- 
modynamic sense an effective trap depth to an ion cloud 
in a multipole ion trap. For this purpose we measure the 
rate of evaporation of trapped ions as a function of their 
translational temperature controlled by a bath of helium 
buffer gas. We use Cl~ anions with a high electron affin- 
ity, because one can safely ignore losses due to parasitic 
chemical reactions with the background gas. The mea- 
sured evaporation rate is found to be proportional to a 
Boltzmann factor with the activation energy given by an 
effective trap depth, i.e. the minimal kinetic energy ions 
need to escape from the trap. We study this effective trap 
depth as a function of the applied rf amplitude and find 
two distinct regimes. Based on numerical calculations, 
the steeply rising trap depth for low rf amplitudes is ex- 
plained by an entirely adiabatic motion of the trapped 
ions, whereas the regime of almost constant trap depth 
for high rf amplitude shows the appearance of regions of 
unstable ion motion inside the trap. An analytical model 
allows us to predict the parameters for a maximum trap- 
ping potential and trapping volume. 

The 22pole ion trap used in our experiment ap- 
proximates a multipole rf field of order n = 11. Storage 
is achieved by 22 stainless steel rods (1 mm diameter) 
forming a 40 mm long cylindrical cage (inscribed diame- 
ter 2ro = 10 mm, see schematic view in Fig.[T]). The rods 
are alternatingly connected to the two ports of an rf oscil- 
lator (w = 27r X 4.7 MHz) providing a cylindical effective 
potential in the radial direction. Along the axis ions are 
confined by small dc voltages (3-10 V) applied to cylin- 
drical entrance and exit electrodes. We manipulate the 
temperature of the ion ensemble via heating or cooling 
the trap and its housing, into which we apply helium at 
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FIG. 1: Measured loss rate of Cl^ anions out of the 22pole 
trap (see schematic view) versus temperature at an rf ampU- 
tude of Vb = 33 V. The sohd hne is a fit to Eq. ((TJ. A similar 
result is obtained with a ID numerical model (see inset) in the 
same temperature range at an rf amplitude of Vb = 10.5 V. 



a well defined density of typically 2 x lO^^cnL""^. Ther- 
malization of the ions hence occurs on the timescale of 
100 lis, assuming a Langevin-limited collision rate. The 
ion source employs a pulsed supersonic expansion of ar- 
gon with a small admixture of CCI4 bombarded by a 
pulsed 1 keV electron beam. In the created local plasma 
Cl~ ions are formed via dissociative attachment of slow 
electrons to CCI4 and transferred to the trap with a time- 
of-flight Wiley-McLaren mass spectrometer. To optimize 
the trapping efficiency ions are initially trapped at a fixed 
rf amplitude for 200 ms, which is then ramped down lin- 
early within 100 ms to the desired value. Typically 2x10^ 
Cl^ anions are trapped per filling, which results in a 
Spitzer self collision time of about 40 s at 300 K, en- 
suring that ion-ion interactions do not play a role. After 
extraction the ions are mass analyzed in a second time- 
of-flight stage before being detected on a microchannel 
plate. 

We observe that the loss rate of Cl~ anions from the 
trap depends strongly on their temperature as can be 
seen from Fig. [TJ For this measurement we heat the trap 
to 370 K and measure the loss rate while slowly cooling 
the trap to 300 K. Each loss rate is derived by measuring 
the decay of the chlorine anion peak in the mass spec- 
trum when increasing the storage time. It is important to 
note that we do not observe any other peaks in the mass 
spectra for any of the storage times. Actually, chemical 
reactions with any of the possible impurities in the back- 
ground gas are not expected considering the high electron 
affinity of chlorine of 3.6 eV. Destruction of the anions 
through detachment of the excess electron by blackbody 
radiation is excluded on this timescale regarding the high 
photon energies needed. Similarly, detachment by colli- 
sions can be excluded in the temperature regime under 
study. This indicates, that the ions are physically re- 
moved from the trap and that the only possible loss pro- 
cess from the trap is evaporation over the effective energy 



barrier formed by the rf field. This is strongly supported 
by the dependence of the loss rate on the temperature of 
the ion ensemble. Fig. [T] reflects that by increasing its 
temperature we force more ions over the barrier. 

We fit the measured loss rate k{T) to a Boltzmann 
factor 

k{T)^Axe^^, (1) 

with the activation energy for loss of anions Ea and the 
pre-exponential factor A as free parameters. Ea is di- 
rectly obtained in units of the Boltzmann constant ks as 
the slope of the fit in Fig. [T] Ea represents the kinetic 
energy ions need to exceed the effective barrier formed 
by the rf field and therefore provides an elegant way to 
determine the effective trap depth. We have determined 
the effective trap depth of the 22pole ion trap as a func- 
tion of the applied rf amplitude Vq. Experimental ac- 
curacies are derived from the Boltzmann fits. As can 
be seen from Fig. [5] the effective trap depth U{Vo) rises 
steeply for low rf amplitudes Vq and reaches a maximum 
of 0.65 eV at around 12 V. For higher amplitudes U{Vo) 
decreases slightly before it eventually levels off at 0.5 eV 
for amplitudes Vq > 40V. This shows that trapping is 
characterized by two distinct regimes, one with a fast in- 
creasing trap depth and one where the depth is roughly 
constant. In addition a small dip in the measured effec- 
tive trap depth is observed at an rf amplitude of 30 V. 

Additional evidence for the two different evaporation 
regimes stems from the pre-exponential factor A in equa- 
tion For large rf amplitudes A is found to be of 
the order of 10^/s, which fits best to the frequency of 
the field, and is independent of the helium density up 
to 2 X 10^^ cm~'^. For very small rf amplitudes A is two 
to three orders of magnitude smaller and compares well 
with the collision frequency of the anions with the helium 
buffer gas. 

To get a microscopic insight into the involved loss pro- 
cesses, we have set up a simple numerical model; details 
will be published elsewhere. Cl~ ions are propagated 
in a one-dimensional oscillating electric multipole field 
of order n ^ E{r,t) = ^ n |r|""^cos(wt + $) with 
—To ^ r < +ro by solving the equation of motion nu- 
merically. We record the loss of ions starting at r = 
with random phase $ with an initial velocity drawn from 
a Maxwellian distribution. The inset in FiglT] shows the 
temperature dependence of the loss rate. In analogy to 
the experiment we determine the height of the effective 
barrier as a function of the rf amplitude. As seen in Fig. 
^jp the numerically determined effective barrier height 
rises steeply for low rf amplitudes, reaches a maximum 
and then slightly decreases for high amplitudes. This be- 
havior is consistent with the corresponding experimental 
data in Fig. [2^. In detail, however, there are distinct 
differences in the maximum barrier height and its cor- 
responding rf amplitude. Nevertheless, the simple ID 
model provides valuable insight. The initial velocities 
of those ions lost during propagation demonstrate that 
exclusively hot ions in the Boltzmann tail of the distri- 
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FIG. 2; Experimentally determined efTective trap depth (a) 
and numerically obtained efTective barrier height (b) as a func- 
tion of the rf amplitude. Each point is derived from a fit sim- 
ilar to Fig. [T] Inserted lines stem from the analytical model 
for different parameters of the maximum adiabaticity parame- 
ter rjmax- The two inserted phase space trajectories represent 
ions with a starting velocity just below the threshold for direct 
escape. At low rf amplitudes the ion returns from the barrier 
with the initial velocity. Strong transfer of energy from the 
field to the ion motion is seen at high rf amplitudes. 



bution are lost, which backs our analysis of evaporation- 
driven ion loss. The calculated ion trajectories in phase 
space are characterized by a uniform velocity in the re- 
gion of low field at smaller radii evolving into a fast os- 
cillation of the velocity in the region of high field near 
the turning point. For a small rf amplitude Vq where in 
Figl2b the effective barrier height is found to be steeply 
rising, no transfer of energy from the field is observed. 
Either the ions are reflected from the flapping potential, 
preserving their kinetic energy, or their initial kinetic en- 
ergy exceeds the effective barrier and ions are lost on first 
approach. The inset in Figl^b for an rf amplitude of Vq 
= lOV shows such a trajectory for an ion starting with 
a velocity just below the limit for getting lost. For high 
rf amplitudes, however, ions with an initial velocity just 
below the limit for getting lost on the first approach are 
found to rapidly aquire energy from the field. This can 
be seen from the trajectory shown in Fig[2l3 for an rf am- 
plitude of Vq = 110 V, where the ion velocity on return 
is more than twice the initial value. The ion will not 
survive the next approach to the barrier. 



By comparison with the numerical simulation we are 
able to interpret the abrupt break-off in the experimen- 
tally observed effective trap depth shown in Fig[2^ and 
the accompanied change in the pre-factor A to a value 
of the order of the rf frequency. We identify it with the 
advent of non-adiabatic ion motion. Ions are lost from 
the trap by rapidly acquiring translational energy from 
the rf field within a few rf cycles, once they happen to be 
located in a trap region where the rf amplitude is strong 
enough to disturb adiabatic motion. 

For an analytic description of the trap depth we intro- 
duce the trapping volume as the volume in which stable 
trajectories without transfer of energy from the field are 
possible. The numerical simulation shows in agreement 
with our measurement that the trapping volume is re- 
duced for increasing rf amplitude and that the reduction 
prevails over the increase in the repelling force by the 
rf field. We parameterize the trapping volume based on 
the adiabaticity parameter r/ defined in ^271 , which scales 
with the (n — 2)-power of the radial position of the ion. 
We consider the trapping volume to be bound by a crit- 
ical radius Vcrit where rj = rjmax if this radius is smaller 
than the geometrical radius rg of the trap, otherwise by 
tq. We then identify the time-independent effective po- 
tential at the edge of the trapping volume with the eflFec- 
tive trap depth. This results in an effective trap depth 
of an ideal multipole trap of 

= (2) 
o e 

where = Mot {ro, T'crit} is the radius of the trapping 
volume, Vcrit the critical radius 

V n-l qVoJ 

and e = 1 / (2n^) mw^rp the characteristic energy. For 
low rf amplitudes Vq the adiabaticity parameter rj does 
not exceed rjmax in the entire trap. Here the trapping 
volume equals the geometrical volume of the trap and 
the effective trap depth is given by the effective potential 
at radius vq. In this regime, the effective trap depth 
grows quadratically with the applied rf amplitude Vq. 
The maximal trap depth is obtained, when rj reaches the 
largest allowed value rjmax just at the geometrical edge 
of the trap. For higher amplitudes, the region of non- 
adiabatic ion motion begins to penetrate into the trap. 
The trapping volume is now bound by the condition 77 = 
Vmax and the effective trap depth is given by the effective 
potential at this reduced radius, which scales as Vjj~^^^. 

The lines in Fig[2^ represent the effective trap depth 
U{Vo) according to the described analytical model for 
Vmax between 0.3 and 0.4. With these values the experi- 
mental data are very well described. Also the numerical 
result in FigO) is described well by Eq. however 
with a higher value of the maximal adiabaticity param- 
eter r]max of between 0.5 and 0.6. This difference might 
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FIG. 3; Effective trap depths for multipoles of different order 
n as derived from the analytical model assuming a maximum 
adiabaticity parameter rjmax = 0.4 (calculation for CP ions, 
ro = 5mm, co = 2-k x 5 MHz). 

be explained by weak fringe fields disturbing the per- 
fect multipole configuration which could originate from 
charge up of insulators and residual gas deposited onto 
the electrodes. This is supported by the experience that 
multipole ion traps provide longer lifetimes immediately 
after cooling down as compared to many days of opera- 
tion at low temperatures. The dip in the effective trap 
depth observed at amplitudes around 25 V is currently 
unexplained and subject to further work. 

The analytical and numerical models can be readily 
applied to multipoles of different order n. In Figl3] the 
effective trap depth for Cl~ ions is plotted as a function 
of the rf amplitude for various multipoles of the same in- 
scribed radius = 5 mm operated at a fixed frequency 
u = 2TT X 5 MHz assuming a maximum adiabaticity pa- 
rameter rjmax = 0.4. The larger field-free region for in- 
creasing multipole order n is payed by a smaller maximal 
effective trap depth. Also, lower order multipole traps 
can be operated in an all-adiabatic mode up to higher 
rf amplitudes. For a quadrupole trap {n = 2) the adia- 



baticity parameter r] becomes independent of the radius, 
but is still a function of the rf amplitude Vq. Note that 
this results in a distinct difference to higher order mul- 
tipoles: The trap depth scales quadratically with the rf 
amplitude Vq until the maximal stability parameter is 
reached; at this point there is a transition from adiabatic 
to nonadiabatic ion motion in the entire trap and the 
trap depth vanishes completely. 

In conclusion, we demonstrate evaporation of an- 
ions out of a multipole rf ion trap. Analysis of the 
temperature-dependent evaporation rate allows us to ex- 
tract the effective trap depth. This shows that Boltz- 
mann statistics are not only applicable to the transla- 
tional and rotational degrees of freedom of the bulk of 
trapped ions [11], but also to the high-energy tail of the 
Boltzmann distribution. Furthermore, we observe the 
transition from an all-adiabatic trapping to energy trans- 
fer from the field to the ion motion. Based on a numerical 
calculation we introduce the concept of a trapping vol- 
ume, and by its parameterization we obtain an analytic 
expression for the effective trap depth and its scaling with 
the rf amplitude, which is directly applicable to multipole 
fields of arbitrary oder n. 

The present results suggest possible applications of 
trap losses as a probing scheme for inelastic collision pro- 
cesses of t rapp ed ions, similar to schemes used for neutral 
atoms (28l.l29|. Maximizing the trapping volume of multi- 
pole traps based on our findings should allow one to trap 
an optimum number of ions at a given space charge inter- 
action, which is particularly interesting for the loading of 
shallow surface traps Furthermore, ion-ion induced 
evaporation in a collision free environment may be inves- 
tigated, for which we have first experimental evidence at 
larger ion ensembles. 
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